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Introduction 

This part of our work focused on understanding the effects of arbitrary transverse and 
longitudinal relaxation rates on the susceptibilities of coherently driven three-level systems. 
The approximation of a single relaxation rate often made in previous work is strongly invali- 
dated by the variation in the spontaneous emission lifetime between various atomic level pairs 
in systems such as cesium. 

It is of great importance to the problem of nonlinear infrared generation to determine the 
dependence both the real and imaginary susceptibility on relaxation rates. The imaginary sus- 
ceptibility on the pump transition determines the absorption of pump photons and the ima- 
ginary susceptibility on the laser transition determines the spectral dependence of the gain. 
This is of particular importance for pure Raman emission (i.e. absorption at linecenter of the 
gain transition) as it determines the tunability characteristics we are aiming to predict. The 
real susceptibility is important when cavities are used at the signal field as this will determine 
the loaded resonance of the Raman oscillator. We will show that in some cases which result 
from having different relaxation rates that mode splitting may result allowing more than one 
frequency to have the same Raman wavelength, possibly resulting in a temporal instability. 
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Study of the Relaxation Rates in Coherently Driven Systems 

The work on elucidating the rate dependence of the susceptibility is based on new solu- 
tions for the three-level density matrix driven by two classical fields of arbitrary strength. 

The density matrix formalism we use is fully equivalent to the Schrodinger wave func- 
tion but is more applicable to the statistical case, in which the wave function is not known 
exactly. In this case the elements of the density matrix are taken to be ensemble averages of 
possible configurations of the system. As a result every wave function can be expressed as a 
unique density matrix, but not all density matrices can be expressed as wave functions, 
reflecting the fact that the wave function contains only quantum uncertainties where as the 
density matrix contains statistical as well as quantum uncertainties. 

For a given wave function 

W> = ZCJn> (1) 

the density matrix is written as 

Pij = cTc; (2) 

where the bar denotes the ensemble average. The time dependence of the elements is given 

by 

Pij = (3) 

where H is the total Hamiltonian of the system used to describe the three level atoms. 

The Hamiltonian in the dipole approximation is given by: 

H = H 0 - frt (4) 

where H 0 is the unperturbed atomic Hamiltonian and -p/S is the atom field Hamiltonian in the 

dipole approximation. We will assume that the dipole matrix element p. 12 is zero correspond- 
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ing to levels |1> and |2> being S-states. The pump transition connects the levels |1> and |2> 
and the Raman emission is near the frequency of the |2>— »|3> transition. 

The transverse and parallel relaxation times per molecule level pair are given by x and T, 
respectively. The levels are assumed to have equilibrium occupation ratios fjj (ratio of i* to j 01 
equilibrium populations). The energy level diagram with the associated rates are shown in 
Figure 1. 

The time evolution of the density matrix elements for the three level system are therefore 
described by 


Pll - T"(P-3lPl3 ~ M^13p3l) ‘^(t) - ~( f 3lP]] - P33)“^7*( f 2lPll ~ P 22 ) 


( 5 ) 


P22 - ~ 7 ~(l-*- 32 p 23 “ P 23 P 32) '^(0 ~ q~( 1 32 p 22 “ P33) - "^~(P22 ~ f 2lPll) 
h Tj T 3 

P33 = |(M-3lPl3 ~ P13P31) + (P23P32 - M-32P23) - — (P33 — f32P22) _ '^ _ (P33 - f3lPll) 

Pl3 = ^ [(P33 - Pll)M-13 ~ + *(©31 + iT2 J )Pi3 

P32 = |(P33 - P 22)^32 - PoiPttj '^(0 ~ i(®32 _ ixr 1 )P32 


pl2 “ h L 


P-13P32 ~ M-32P13 


j-S(t) + i(co 2I 


+ *^3 ')Pl2 


( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 

( 10 ) 


where C0jj = (E, - Ej)/h, E(t) = E p cos(co p t + \p p ) + E s cos(co s t + \j/ s ) and p 31 , p 23 and p 12 are complex 
conjugates of pi 3 , p 32 and p 12 respectively. 


Steady State Solutions 


To obtain the steady state solutions to the density matrix equation we first take out the 
explicit oscillations of the off diagonal elements 
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P \3 = Ae ia3pt p 31 =^‘ iQ3 ‘ l p 12 = De i(Up ~ Ws)l (11) 

and then neglect all terms which oscillate faster than co p or co s in the off diagonal element 

equations and all non DC terms in the diagonal element equations. This is the rotating wave 
approximation which therefore disallows the possibility of Bloch-Seigert shifts in the reso- 
nances. If we then set all time derivatives equal to zero and define the following: 


P s = 


M-23 E s6 

2h 




Pp = 


P-I3 E p e 

2h 


‘V. 


( 12 a) 

( 12 b) 


L s = ~ ffl 32 + i/Xj = 5 S + i/t] (13a) 

L p = co p - co 31 - i/t 2 = 5 p - i Jx 2 (13b) 

L sp = co s - co p + co 2 i + i/t 3 = 5 S - 5 p + i/x 3 (13c) 

and A 32 = p 33 - P 22 A 13 = P 11 _ P 33 the equations become 

L s \=P s *A 32 -P;d (14a) 

LpA = -P p A 13 - (3 S D (14b) 

L sp D = -p p X + p*A (14c) 

and 

0 = 2ImP s \ - 4ImP p *A + k^A,, - A, 0 ,) - k 2 (A 32 - A 3 ° 2 ) (15) 

0 = 2Imp* A - 4Imp s X + k 4 (A 32 - A 3 ° 2 ) - k 3 (A 13 - A [ 3 ) (16) 


The constants k 1# k 2 , k 3 , k 4 , A^ and A 32 are defined in terms of the T and f ; j by the following 


relations: 
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We can now eliminate D from the first three equations and solve for X and A 


where 


-14 




a 32 + IP p ! R 


M3 


+ 

L L p L i 


32 



(17) 

(18) 


R = 


Lsl^p 


LsiP s i 2 -L p U3 p | 2 + L s LpL sp 
In terms of two real quantities, R[ and R 2 , we have that: 


R = Rj - iR 2 = 


(5 2 + VxS^ + VXA-iB) 

A 2 + B 2 


_ w P 


A = 5 s (8 p + t 2 2 )(5 2 + Xi 2 - !p p | 2 ) - 8 p (5 s 2 + Xj 2 )(5 p 2 + x 2 2 -. |p s | 2 ) 
B = x 3 -1 (5 2 + xf 2 )(5 p + x 2 2 ) + lP p lV(8p 2 + x 2 2 ) + IPslViS 2 + xf 2 ) 


This then yields 


Imp s 7, = 


-!Psl : 

|L S | 2 


iPpt 2 


V + 


|L S |‘ 


a 3 


A32 - 


IP P I 2 IP S 1 : 
I L s 1 2 1 i 


' a 2 Ai 3 


(19) 

(20) 

(21) 

( 22 ) 

(23) 


and 
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Imp p *A = 


- 1 Ps r 
ILpI 2 


,-i iPpi 2 

2 |L p | 2<Xl 


Ai 3 — 


IPpI 2 IPsI 


14414 


2 i t ,2^32 


where 


«i = 2t 2 " 1 5 p R 1 - (8p - x 2 2 )R 2 
02 = (Sptf 1 - 8 s t 2 “ 1 )R 1 + (xf’xj 1 + 8 p 8 s )R 2 
a 3 = 2xf 1 5 s R 1 + (8 S 2 - xr 2 )R 2 


Substituting these into the last two density matrix equations then gives 


44°3 ~ 4^32 - ~ 


2 IPs!' 


141' 


+ < 


4iPp| 2 


14 


V - IM : 


-rj 1 + 2|p_| 2 


Oj «2 

T 


-o 2 a 3 


ILpI 2 I4I 2 


+ k-i 


4 


32 


141 2 | 4 |' 




+ k. 




13 


4^32 “ 4^13 - + 


4|Ps 


141' 


xf l + IPpi' 


-a 2 a 3 


2 i L p 14 P 


2 1 P p 1 2 


14 


x 2 1 - IP. I 3 


Oi 2a 2 

+ 


141 141' 


+ ki 




13 


Then by making the following substitutions and definitions 


+ 4 


4 


32 


4 = 4+ _ 2 


4 1 P p 1 2 
5 P 2 + xt 


V “ IPs!' 


Oi 


+ 


a 2 


L 5 p +- x 2 z 2(8 S Z + xf 2 ) 


r 2 = ^fr- + 2, P ^ 


-a 2 


8p + x 2 2 2(8 S 2 + xf 2 ) 


+ 


<*3 


■ + k-j 


rj . *!ediL.- l(y » 

4 +^2 2 I 


a l 


2a? 


5p +x 2 z 5 2 + xf 2 


■ + fc> 


(24) 

(25a) 

(25b) 

(25c) 

(27) 

(28) 

(29a) 

(29b) 


(29c) 
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r. = k.+ 4|Psl :? |rr l + !P p i 2 


5 S + Ti ■' 


-02 03 


2 ( 5 p + t 2 2 ) 5 2 + x , 2 


The equations become 


with the solutions 


A13 - 


A 32 - 


kiA°3 - k 2 A° 2 - TjAn - T 2 A 32 
^4^32 — ^4^l°3 = -r 3 A,3 + r 4 A 32 

r 4 (kiAf3 — k 2 A 32 ) + r 2 (k 4 A 3 ° 2 — k 3 AP 3 ) 

r,r 4 - r 2 r 3 

ri(k 4 A 3 ° 2 - k 3 Ai°3) + r 3 (k 1 A 1 Q 3 - k 3 A3 2 ) 

r,r 4 - r 2 r 3 


(31a) 

(31b) 


The polarization of the system can now be determined from the expectation value of the 
dipole operators by: 


<P> = Tr(pp) 

= 2Re(p 31 p 13 ) + 2Re(H 32 p 23 ) 
which, after substituting for p 13 and p 32 becomes 


P = 2Re 


'^ 13 T“ i A 13 - I Ps I 2r 


2hU 


+ 2Re 


_ I M-23 I ‘ 


2hL c 


a 32 + I P p I Z R* 


A 13 A 32 

l„ u 


A 13 A 32 

L* L* 
L l p 


rE s e 


i(w p t+\|/p) 


i(io s t+Vj) 


By comparing this equation with 


P = Re 


XpE/'^V] + Re (^W 


the complex susceptibilities are found to be given by: 


x P =d sr" A ' 3 “ IPi|2R 


A 13 A 32 

L L P L s 


( 32 ) 


(33) 


(34) 


(35) 


(29d) 

(30a) 

(30b) 


and 



Progress Report/NASA 


- 8 - 


NAG 5-526 


+ I M-23 I ' 


A32+ lPpl R TT + T^ 


In their complex forms 


X P = X P ' + iX P 


Xs = Xs + iXs 


~ I M'H 1 2 f L-i IPsl 2 «i 
Xr h(8 2 -Kf 2 ) 1 13 [ 2 (Sp + X 2 - 2 ) 

Ite l : 2 [a L-i , !Pp |2 °^ 

h(5 s 2 +xf 2 ) I 32 [ 1 (8 S 2 + xf 2 ) 


X P ' = 


-ip- 131 2 ( L _ 1 P s 1 2 yi 

h(8 2 + xf 2 ) 1 13 [ p (5p + xf 2 


1^32 1 2 r l , iPpi^ 

* s WX 2 , „-- 2 n r 32 5s + /s 2 , _- 2 i 


h(8 s 2 + * 1 2 ) 


IPsl 2(x l 1 , . IPsl 2 «2 I 

Sp + t 2 “ 2 ) 32 (8 2 + Xf 2 ) J 

I Pp I 2(x 3 1 A I P p 1 2 Ct 2 ] 

:s 2 + xf 2 ) 13 (Sp + t 2 2 ) J 

1PsI 2 Yi 1 IPsl 2 Y2 

(Sp 2 + T 2 2 ) J 32 (8 2 + xf 2 ) 

lPpi 2 Y3 +A IPpl 2 Y2 
(S s 2 + xf 2 ) _ 13 (8 2 + xf 2 ) 




(36) 

(37) 

(38) 

(39) 

(40) 


Yi = (S 2 - t 2 “ 2 )Ri + 2x 2 - , 8 p R 2 (41a) 

Y 2 = (5 P 5 S + xf’xf *)Rj + (8 s x 2 _i - 5 p xf 1 )R 2 (41b) 

Y 3 = (S 2 - Tf 2 )R! - 2xf>8 s R 2 (41c) 

These susceptibilities are the first results which allow for completely arbitrary relaxation 

and dephasing rates as well as arbitrarily strong pump and signal fields. 

The systematic study of the relaxation rate dependence was performed by setting two of 
the pairs of rates (longitudinal, T and transverse x) equal and allowing the other pair of T ; and 
X; to have a different ratio. This amounts to eliminating the unity ratio of transverse to longi- 
tudinal rates between one pair of levels at one time. The results were astonishing and indicate 
that the. single ratio limit is a pathological case of no great value for most atomic systems. 

The next section will give results due to the variation of the ratio of longitudinal to 
dephasing times. The curves will be labeled in terms of 8 P and 8 S ; the pump and signal field 
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detunings in terms of T" 1 units, and p P and p s the pump and signal field strengths in terms of 
saturation field values. For example, a p P = 1 indicates a value of pump intensity equal to the 
saturation intensity based on a two-level model. 



Xp fi /^ 2 e 0 Xpfi/)i2 
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Pump Susceptibilities %p", % P ' 


8 S — — 5 |3p — 2 Ps — 2 
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Mode Splitting and Relaxation Rates 

It has been realized since the discovery of pulsations in inhomogeneously broadened Xe 
lasers that complicated lineshapes (hole burning in the case of the Xe laser) could lead to con- 
ditions where a single cavity mode could be consistent with several frequencies. This is stat- 
ing the fact that owing to complicated frequency dependent dispersion, several frequencies can 
have the same wavelength. In our study of the signal susceptibility dependence on the T/x 
ratio we discovered that the possibility of mode splitting due to nonunity ratios exists and is 
likely to result in transient pulsations in the Cs system if a cavity were employed. As an 
example we give the graphical solution of the Fabry-Perot Cavity resonance equation 


1 + 


for two relaxation rate ratios for the case of stron 


L = m;t 

g off-resonance pumping of a three-level 


system (5 P = -5, P s = 0, P P = -2): 
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The intersection of the straight line with the X / (6g) curve locates the frequencies where a cav- 
ity tuned to the center of the 1 3>- 1 2> transition would initially begin to oscillate. The curve 
shows that these frequencies (labeled co, co l5 CO 2 and co 3 ) can satisfy the oscillation condition 
when T 2 = lOt while only one can oscillate when T = x. 

Conclusions 

This phase of our work has clearly shown that all previous models for Raman emission 
in coherently driven three-level systems are in most cases invalid for (1) predicting spectral 
signatures and (2) predicting Raman oscillation frequencies. The various cases in the text 
show that the ratio of the transverse (dephasing) to longitudinal (population) rates on any pair 
of levels in the three-level system can strongly affect the susceptibilities. These situations are 
in fact the rule rather than the exception when atomic species such as Cs, K or Na are con- 
cerned. 

In addition we have shown that when strong off resonant pumping such as would be 
required for tunable JR generation, a larger dephasing rate results in strong modulation of the 
real susceptibility. Using a graphical solution for the threshold resonance of a Fabry-Perot, 
we have shown that such a Raman oscillator could begin oscillating on several frequencies. 
This would most likely lead to temporally varying output even for smooth injected pulses. 

The work we are currently pursuing involves the inclusion of the doublet effects present 
in the Cs upper state (|3>) on the x'' and %' calculations. We will shortly produce predictions 
of the time averaged IR output as a function of the power and frequency of the pump for 
comparison of the experiments at Goddard Space Flight Center. 



